arXiv: 1508.0687lv2 [math.NA] 12Jun2016 


Discrete Green functions of the SDFEM on 
Shishkin triangular meshes 

Jin Zhang*^ 

June 14, 2016 


Abstract 

We present bounds of discrete Green functions in the energy norm 
for the standard (or modified) streamline diffusion finite element method 
(SDFEM) on Shishkin triangular meshes. 


1 Problem 

We consider the singularly perturbed boundary value problem 

—eAu + b ■ Vu + cu = f in O = (0,1)^, ... ... 

u = 0 on dfi, ' 

where e « |6| is a small positive parameter, b = (bi,b 2 )'^ is a constant vector 
with 6 i > 0 ,62 > 0 and c > 0 is constant. It is also assumed that / is sufficiently 
smooth. The solution of (El) typically has two exponential layers of width 
0 (eln(l/e)) at the sides x = 1 and y = 1 of fl. 


2 The SDFEM on Shishkin meshes 


2.1 Shishkin meshes 

When discretizing dm, we use Shishkin meshes, which are piecewise uniform. 
See [SJ [3 [5] for a detailed discussion of their properties and applications. 

First, we define two mesh transition parameters, which are to be used to 
specify the mesh changes from coarse to fine in x— and y—direction. 


A 3 , := min 




1 6 

and Xy := minip—lnN 

2 P2 
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Figure 1: Dissection of D and triangu- Figure 2: K^j and Kfj 

lation Tat- 


where /3i and P 2 are defined as in [31 Assumption 2.1]. For technical reasons, we 
set p = 2.5 in our analysis which is the same as ones in |S] and [5]. The domain 
is dissected into four parts as D = u u u (see Figure [1]), where 


ns := [0,1 - AJ X [0, 1-Xy], n^ := [1 - A„ 1] x [0,1 - Xy], 

ny := [0, 1 - AJ X [1 - Xy, 1] , n,y -.= [l - A„ l] X [l - Xy, l] . 


Assumption 1. We assume that e ^ N~^, as is generally the case in practice. 
Furthermore we assume that Xx = In A and Xy = pefF^^ In A as otherwise 

N~^ is exponentially small compared with e. 


Next, we define the set of mesh points 


I 2i(\-Xx)IN 
1 l-2(A-f)A,/A 


{{xi,yj) e n : i, j = 0, ■ ■ ■ , A} 

for i = 0, • • • , A/2, 
for i = A/2 + 1, • • • , A 


and 


I 2j(l - Xy)/N for J = 0, • • • , A/2, 

1 1 - 2(A - j)Xy/N for j = A/2 + 1, • • • , A. 


By drawing lines through these mesh points parallel to the i-axis and j/-axis 
the domain Q is partitioned into rectangles. Each rectangle is divided into 
two triangles by drawing the diagonal which runs from {xi,yj.^.l) to 
This yields a triangulation of D denoted by 7Ar(see Fig. [T]). The mesh sizes 
hx,i ■= Xi+i - Xi and hy^ := pj+i - yj satisfy 


hx.i 


Hx := 


l-Xx 

A/2 


hx • — 


A/2 


for f = 0, • • • , A/2 - 1, 
for i = A/2, • • • , A - 1 
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and 


for j = 0, ••• , Af/2 -1, 
forj = iV/2, ... ,iV-l. 

The mesh sizes and hyj satisfy 

N-^ ^ H^,Hy ^ 2N-^ and CieTVMnTV ^ h^,hy ^ C'2eiV“Mn7V. 

For convenience, we shall use some notations: Klj for the mesh triangle 
with vertices {xi,yj), {xi+i,yj) and (xi^yj+i); Kfj for the mesh triangle with 
vertices (xi^yj+i), {xi+i,yj) and {xi+i,yj+i) (see Fig. [2|); t or K for a generic 
mesh triangle. 

2.2 The streamline diffusion finite element method 

On the above Shishkin meshes we define a linear finite element space 

:= e Cifl) : v^\gn = 0 and v^\k e PiiK), VK e Tn}. 

Now we are in a position to state the standard SDFEM for dm which reads: 
f Find e such that for all e , 

I asoiu^,v^) = {f,v^) + Yj U,^Kb-Vv^)K, 

I KczQ 


where 


asoiu^, v^) =e(Vu^, Vv^) + (b ■ Vm^, v^) + {cu^,v^) 

+ ^ {—eAu^ + b ■ + cu^, 6Kb ■ Vv^)k 


Note that A{u^\k) = 0 for u^\k e Pi{K). Following usual practice [7], the 
parameter 5k '■= 6\k is defined as follows 


6 := 6{x,y) 


(C*N-^ ifix,y)Bn,, 
\ 0 otherwise, 


( 2 . 2 ) 


where C* is referred to [7l Lemma 3.25]. 

We set 

b-.= ■\Jbl + bl, r] /b and 


for any vector of unit length. By an easy calculation one shows that 


(Vw, Vu) = {wp.vp) + {Wy,Vy). 
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We rewrite (EH) as 


+ e{u^,v^) + {bu’^ + u^,v^)+ ^ {bu^ + cu^,5Kbv^)K 
= (/w'^)+ E {f,SKbv^)K. 


KtzQ. 


Ktzn 


For technical reasons in the later analysis, we increase the crosswind diffu- 
sion(see [TJI31E]) by replacing e{u^,v^) by i{u^,v^) where e ^ e ^ N~^ and 
e is constant on each of subdomains including ns,fl\fls- For convenience, we 
denote e|n,by fs- We will consider two cases: (1) e = e\ (2) e defined as in [3l 
pg 463], that is 


e: = 


if a; G 

if a; G Q\'[ls 


where e := max (e, . 

We now state a streamline diffusion method with artificial crosswind diffusion 
(ACD), also called a modified SDFEM: 


Find G such that for all G 
aMSoiu^, v^) = (/, + Sbvp), 


with 


aMSD{u ,v^) :=e{up,vp) + e{u^ ,v^) + {bup +u^ ,v‘^) 


+ Xj ,5KbVp)K- 


(2.3) 


Kczd 


Finally, we define a special energy norm associated with UMSoi’, •)• 

Ik^llMsn •= + Xj yv^eV^. 


KczQ 


For brevity, we often write e(u^, v^)+ ^ {bu^^Sx bvp)K 3.S {e+b‘^S){up ,Vp). 

KczQ 


3 The discrete Green function 


Let a:* be a mesh node in ft. The discrete Green’s function G G associated 
with a;* is defined by 

aMsniv^.G) = {x*) g . (3.1) 


The weight function to is defined by 


uj{x) ■= g 


{x — X*) • / 9 ^ f ^ X*) ■ V 
crp 


9\- 


(x — X*) ■ r] 
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where 

2 

g{r) = for r e (- 00 , 00 ). 

Now, we are to derive a global estimate on G in the weighted energy norm 

111^111^ :=e||a; + e||a; +-||(a;“^)^'^^Gp (3.2) 

+ cWuj-^I^GW^ + ^ 95k\\uj-^'^Gp\\. 

K 

From (|2.3p and (|3.2p . we have 

IIIGIII^ =aMSD{uj-^G, G) - e{{u;-^)pG, Gp) - i{{u;-\G, G,) (3.3) 

- 2(6(o.- 1)^G + cw-^G, Sk bGp)K- 

K 

Considering (13.11) we also have 

aMSD{uj~^G, G) = aMSD{E, G) + aMSD((w“^G)^, G) 

= aMSD{E,G) + ioj-^G){x*) 

where E := uj-^G - {uj-^GY . 

Lemma 3.1. If 

1 yo 

ap ^ k{e + Sm) and ct,, ^ ki^ , 

where 6 m ■= max^j^o (5 and em '■= maxajgoe, then for k > 1 sufficiently large 
and independent of N and e, we have 

aMSD{u^-^G,G)^^\\\G\\\l 

Proof. Holder inequalities, Cauchy inequalities and [3l Lemma 4.1] give 
{e + b^S)\{{u;-YpG,Gp)\ 

+ ^G64/V-'/"||(w-i)^/'G|k • bdY^\Y-^^^Gp\\K 

K 

^^'^GpY + Ceap^Yu} + -'^h'^5K\a} ^^‘^GpYx 

^ * K 

+ GY,5KO-pY{^-YfGf 

K 

^^'^GpY + - 'Y^h'^5K\oj ^/‘^GpYk + G{e P 5)ap^\{Lo Y^p'^G^. 

^ ^ K 


(3.4) 


(3.5) 
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Similarly, we have 

e|((a;“^)^G,G^)| ^ Ge^/^cr”^ • ||a;“^/^G|| • e^/^||a;“^/^G^|| 

i£||u;-V2G^||2 + ]^CEa-2 d^-V^Gf 

and 

^(cW-1G,(5k hGp)K < c 1/2||^-1/2^||^ . ^,4/2||^-l/2Gd|^ 

K K 

^ ic||a;“^/^Gp + '^c6k ■ b‘^6K\\i^~^^'^Gj3\\j^ 

K 

^ic||..-V^Gp + i2&^<5;,||u;-V2G,||^. 

^ ^ /c 

If (13.511 holds true and k is taken sufficiently large and independent of N and 
e, then we have 

(e + 6^5)|((w ^); 3 G, G; 3 )| + £:|((a; ^)^G, G^)|+^ |(cw ^G,6k bGp)K\ ^ jWlGlW^. 
Considering (13.31) . we are done. 

□ 


Lemma 3.2. If 1 ^ ap ^ fe(e + Sm), with k > 0 sufficiently large and indepen¬ 
dent of N and e. Then for each mesh point x* G Iffllxy, we have 

|(w-iG)(a;*)| ^ il||G|||2 + ^ ^ 


16 


[GiVlnlV ifx*enxuny 




where G is independent of N, e and x*. 

Proof. First let x* G Us- Let r* be the unique triangle that has x* as its 
north-east corner. Then 

|(a;-iG)(a:*)| ^ GN\\G\\r* < GfVmax 

Calculating {u}~^)p^{x) explicitly, we see that 

(u;“^)^^(a:) ^ Gcr^ VajGr*. 

Thus the arithmetic-geometric mean inequality gives 


\(uj-^G){x*)\^GN^ap + - 


2 

id ■ 


Next, let G Qj.. (The case g 17y is similar.) Write a;* = (xj, yj). Then 


|(i^-'G)(a,*)| = |G(x*)l - 


r 


Gx{t,yj)dt 


^ GH. 


n Vj + i 

\Gxit,y)\ dydt 
h 


^ GN{e\nN ■ N-^y^^\\Gx\\a^ ^ GiV^/^ In^/^ iV|||G||| 
1 


^ GA^lniV + 


16 


2 

oj ■ 
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□ 


Thus, we are done. 

Lemma 3.3. If 

where k > 1 sufficiently large and independent of N and e, then 

aMSD{E,G) ^ 

where E = {uj~^G)^ — oj~^G. 

Proof. Holder inequality gives 

\aMSDiE,G)\ + b^Sy^^Wcu^^^EfiW ■ (e + b^Sy^^Wco-^^^GpW 

+ G\\uj^^'^E\\ ■ \\oj-^l‘^Gp\\ + \\uj^l‘^E\\ ■ \\uj-^l‘^G\\. 

^ — |||G|||^, apply (13.61) and take k sufficiently large. 

16 

□ 


Set hr = max{ha;_T-, hy^r}- For Vr e Tjv we have 

sg G(maxa;)F2 . (||(a;“^G);3;3||^ + ||(a;"^G);3^||^ + ||(w“^G),,,,||^) 

T 

and 

\\uj^^'^Ep\\r ^ G(maxa;)^/^ • hr (||(w“^G)/3/3||r + \\{ui~'^G)/3y\\r + ||(w“^G),,,,||,-) , 
Wuj'^^'^EyWr ss G(inaxa;)^/^ • hr (||(a;“^G)/3/3||T + ||(a;“^G);3,,||r + ||(a;“^G),,,,||r) , 

T 

where we have used [Sj Corollary 3.1]. Clearly, the following estimate holds 

||(w ^G)/3,g||r + IKw ^G)/3,,||r + ||(w ^G),,,,||r 

^C{\\{uj-^)ppG\\r + \\{0J-^)pGp\\r) 

+ G(||(a; ^)/3,,G||r + IKw ^)/3G,,||r + ||(a; '^)yG^\r) 

+ G(||(a; ^),,,,G||r + ||(w ^),,G,,||r) . 

Set a^: := (rnino;)”^/^. Lemma 4.1 in [5] yields 

T 

\\{u:-^)ppG\\r ^ Ga^ ■ • \\{u:-^)fG\\r. 

\\{<w-^)pGp\\r ^ Ga^ ■ ■ ||w-F2g^||,, 

||(cc-i)^,G||. < Ga. • ■ \\{w-^)fG\\r. (3-7) 

\\{oj-\Gp\\r ^ Ga^ ■ u-i • \\uj-^^^Gp\\r, 

||(w-i),,,,G||, ^ Ga^ ■ a-^ ■ ||c^-'/'G||,. 
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Similarly, for r c fig we have 


1 (w ^)^Gr;|T ^ CUuj ■ 

j g-^N -Wuj-^/^GWr Hig^N-^ 

|cr“^e7^^^ • Es^^IIw^/^G^IIt- if £s > N-^ 

(3.8) 

and for r cz n\ng 

(w ^)r)Gni\ 

\riiGau:-e ^ • £^/^||a; ^"^^G^It-. 

(3.9) 


For T cz rig, we have 


||(a; ^ C'max(w ^)/3 • ||G^||t-^ Gmax(a; ^)/3 ■ fV • ||G||^ 

T T 

/ \ 1/2 / \ 1/2 
^ G/V (^max(w“^)^J • I|G||t 

^ Ca^ ■ ■ ||(w-1)^/'G||, (3.10) 

where we have used [3J Lemma 4.1 (vii)]. For r cz n\ng, we have 

\\{u;~^)l3Gn\\T < Ga,^ • • e^/^||w“^/^G^||r. (3.11) 


Set 


7 (e) = 


v-i, 


if £s ^ /V ^ 

if £s ^ N~'^ 

From (I3.7I) - (I3.11I) and [3j Lemma 4.1 (iii)], we have 


g^-N 
-1 --1/2 

G^^Es 


(3.12) 


||wi/2G||n. + iV-i 

^CN + ^/3 + ^/3 ^ + 'I'r? ^ + ^r;(£) + ^/3 III ^111 


l|w'/"G||n\o, + iV-i (||u;i/2G^||o\n. + ||w'/"L;,||n\o.) (3.13) 

^CN +Gp^e + cr^ '^CT^^+CT^^£ + cr^|||G|||,^_ q\o^. 

To make sure that (13.61) holds, we should have 

||c^^/'£//3||n, <Gfc-i/2(e + 5^)-i/2|||G|||,,_o., (3.14) 

||w'/"£^,||a. < Gfc-i/2£7i/2|||G||U,a,, (3.15) 

l|c^'/'i^/3||n\n. + l|w'/'i?,llo\n. < Ck-^^^e-^^^\\\G\\\ ui,n\n^, (3.16) 

^ Gk-^^^Sl^^\G\\\^,n., (3.17) 


and 

l|w^/^G||n\n, < Gfc-i/2£i/2 III Gill uj,Q\Q,s ■ 

To ensure that (13.141) holds true, we can set 


Gp ^ kN 

Gr, > kN-^!'^, 


Ut) > uir := 


J /cAi-i /2 


if £^ < Ai-2 
if £^ ^ Ai-2 ■ 


(3.18) 


(3.19) 
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Assume that (13.191) holds true, then (I3.14|) - (I3.17D hold and we have 


^ (3.20) 


Next, we are to obtain sharper bounds for . Following the tech¬ 

niques of (see [21 Lemma 4.4]), we have 





(3.21) 


where x e n\ns, r(a;) e dil satisfies {x — r(a;)) • /3 = 0 and the following 
condition: 


\x — r(a:)| = minja: — y|, where y e dft and {x — y) ■ {3 = 0, 


From (I3.21L we have 


w 




2 

Q,X^^X 


fl fl r fr(a;) 

—Afr 


n 2 


(a;^/^i?)^ds 


dj/da; 


^ C\l {||(u;'/")^F;||^\a. + } 

Similar argument holds for fly. Thus, we have 


(3.22) 


To make sure that (|3.18ll holds, according to ()3.22p the following estimates 
should hold 


ln-2 iV|||G|||5(3.23) 

From (13.231) and (I3.20F we can set 

up ^ kN-^ In N, Ur, > kN-^ In N, Ur, > ke^l^N-^l"^ In^/^ N. (3.24) 

Assume that (I3.19P and (13.241) hold true, substituting (13.201) and p.23l) into 
the right-hand side of (I3.22F we have 

< Ck-^e\\\G\\\l. (3.25) 

Thus, we have the following lemma. 
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Lemma 3.4. Assume that ap and Un satisfy (13.191) and (13.241) . where k > 1 
sufficiently large and independent of N and e, then we have 




\\uJ^^^E0\\n^ + WujE^E^Wn^ ^ Gk-E^N^^^\\\G\\\^,n., 

+ ||u;i/2L;,||o\n, < \n-^ N\\\G\L^n\n^. 

Theorem 1 . Assume that up and Ur/ satisfy (13.51) . (13.191) and (13.241) . where k 
is chosen so that Lemmas roro and \3.A hold. Then for x* G D\Llxy, we have 


||G||msd < V8|||G|||cc; ^ 


GiVaf 

C'lVi/2 ini/2 TV 


if X* G Lis 
if X* e Ltx u fly 


Proof. The readers are referred to [Sj Theorem 4.1] for the estimate ||G|| Msn < 
V8|||G|||(^. Considering (13.3L (13.41) and Lemmas I3.1H3.31 we obtain 

JlllGIII^ ^aMSD(w-'G,G) = aMSD{E,G) + {w-^G){x*) 

1 lll^ll|2 , J GA^V;3 ifa;*Gfls 

-glllt^llL + |(T7TviniV ii x* e Ll^'u Lly ' 


Thus we are done. 


□ 


From (13.5F (I3.19|) and (13.241) we have 

up ^ kN-^nN, Uy ^ kinax{e]^^,N-^^'^,u*,N-HnN,eE^N-E'^\n/'^N} 

where cr^ = i - 1/2 -j/o 9 ■ Note that if ^ 4, then ^InN 

^ [kis'N-^/^ liisl^N-^ 

and 

7V-1/2 > iV-Hn N, 7V-^/2 ^ gi/4Tv-i/2 for TV ^ 4. 

Now, we consider the case of e = e. Clearly, em = is = s. 

• If e ^ N~'^, then u* = kN~^/'^ and we have 

TV-1/2 ^ gl/2^ tV- 1/2 ^ iV-3/4, 

N-E^ > N-^ In N for iV > 4, 

Tv-1/2 ^ £1/4tv-i/ 2 foi/2 TV for AT ^ 4 . 
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• If TV ^ < e ^ TV then a* = ke ^/^TV and we have 

e-l/2^-3/2 ^ gl/2^ ^_i/2^-3/2 ^ ^-3/4^ 

It means that tr^ ^ /c max{£“^/^TV“^/^, TV“^ In TV, e:^/^TV“^/^ In^^^ TV}. Note 
that if TV-2 < £ ^ TV-3/2, 

TV-1/2 > e-i/2jv-3/2 

iV-1/2 > TV-1 ^ for TV > 4, 

^-1/2 ^ ^l/4^-l/2 j^l/2 ^ for iV ^ 4. 

• If TV-3/2 < e ^ TV-1, fjrgo cr* = and we have 

£1/2 > £-1/2TV-3/2, eV2 ^ fV-3/4. 

It means that cr,, > fc max{£i/2, TV-i In TV, £i/^TV- 1/2 jrii/2 JV}. Note that if 
TV-3/2 < ^ 

TV-1/2 > £1/2 

iV-1/2 > AT-i In TV for TV > 4, 

^- 1/2 ^ gi/ 4 ^- 1/2 jj^i/2 ^ for iV > 4. 

These mean that in the case of e = e we can set tr,, ^ kN~^^^. 

^ ^ othe^w^se ^ max (£, TV-3/2), clearly, em = is = e. 

Note that cr* = fc£-i/2TV-3/2 f^^. ^ ^^- 3 / 2 ^ 

• If £ ^ then em = £ = TV-3/2 have 

42 ^ gl /2 ^ ^-3/4^ = tv-3/4. 

It means that ct,, ^ /c max{ei/2, TV-i InTV, £i/^TV-i/2 Ini/^ TV}. Note that 

£1/2 Ini/^ > TV-1 In Af > 4, 

£1/2 lni/2 A^ > £1/4 tV-i/ 2 lni/2 N hr N ^ 4. 

• If TV-3/2 < £ ^ TV-1, fjrgri e = E > and we have 

£' 1/2 > £-1/2TV-3/2, 

£1/2 lni/2 A^ > A^-i In AT for Af > 4 

£1/2 lni/2 N > e^I'^N-^I'^ Ini/^ A^ for A^ > 4 
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These mean that we can take arj ^ N. 

From the above derivations, we have the following remark. 

Remark 3.1. If we set e = e, we can take 

ap = kN-^\n.N, Ur, = kN-^'^. 


If we set £ as in m pg 463J, i.e., 

. ^ if X ells, 

1 e otherwise, 

where i := max (e, we can define 

ap = In N, cr^ = N. 

It seems that the definition of ap in is not proper. 
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